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Abstract 

The e n -method Is employed with the spatial 
amplification theory to compute the onset of tran- 
sition on a swept wing tested In transonic cryo- 
genic flow conditions. Two separate eigenvalue 
formulations are used. One uses the saddle-point 
method and the other assumes that the' amplification 
vector Is normal to the leading edge. Comparisons 
of calculated results with experimental data show 
that both formulations give similar results and 
Indicate that the wall temperature has a rather 
strong effect on the value of the n factor. 

1 . Introduction 

In the absence of leading-edge contamination 
and Gortler Instability, It Is well known that 
transition on swept wings may occur either due to 
streamwlse Instability (related to the properties 
of the streamwlse velocity profile u) or cross- 
flow Instability (related to the properties of the 
cross-flow velocity profile w) . Since the 
u-proflles look like Falkner-Skan profiles, the 
streamwlse Instability is similar to that of a 
two-dimensional flow and leads to turbulence In 
flows with positive pressure gradient. On the 
other hand, the development of a cross-flow pro- 
file Is characterized by an Inflectional Instabil- 
ity that can Induce transition In flows with a 
negative pressure gradient. 

In order to design new aircraft wings, one must 
be able to predict transition on swept wings. The 
most popular method for predicting transition Is 
the e n -method, which was Initially developed for 
two-dimensional flows^**. j^ s method Is based 
on the solution of the Orr-Sonmerfeld equation 
using either temporal or spatial amplification 
theory. In either approach the Integrated ampli- 
fication rates A/A 0 of the unstable frequencies 
are determined, and transition Is computed on the 
assumption that It occurs when the ratio A/A 0 of 
the locally most unstable wave reaches a critical 
value e n , with n between 8 and 10 for a low dis- 
turbance environment. 

In the extension of this method to three- 
dimensional flows, both temporal and spatial 
amplification theories can again be used. In the 
former case, the eigenvalue problem Involves five 
scalars a, $, « r , and R and In the latter case 
It Involves six scalars a r , a^, 3 r , , w and R. 
In both approaches, the solution procedure Is con- 
siderably more difficult than Its counterpart In 
two-dimensional flows because the nontrivial solu- 
tion of the linear stability equations In three- 
dimensional flows provides only two relations 


between the eigenvalues a, p, u and R. To predict 
transition for example, In spatial amplification 
theory, u and R are prescribed, so two new rela- 
tions between a and p are required before the solu- 
tion of the linear stability equations can be 
obtained. 

To date, most problems of the three-dimensional 
transition problems employing linear stability 
theory have been treated by using the temporal 
theory. When R and « r are specified, the values of 
a and are not unique, since they both are func- 
tions of p [a = a(p) and A possible 

solution Is to determine the wave number direction 
* tan -1 (p/a) for which Is maximum and Inte- 
grate ui along the group velocity direction accord- 
ing to Caster's transformation 3 . This "envelope 
method" Is one of the approaches used In the COSAL 
code developed by Malik. * 

In spatial theory, the Introduction of an addi- 
tional scalar In the dispersion relation makes the 
problem more difficult: the amplification rate vec- 
tor A Is no longer a scalar; It Is a function of 
both a* and 0). As a result, a new relation Is 
needed In the eigenvalue formulation. 

In this paper we consider two completely dif- 
ferent eigenvalue formulations and compare their 
predictions with measurements. The first formula- 
tion Is based on the wave packet theory (saddle- 
point method) and the second Is based on the as- 
sumption that the amplification vector A Is normal 
to the wing leading edge. The experimental data 
correspond to measurements obtained at ONERA/CERT 
for a transonic swept wing. The tests have been 
conducted In a cryogenic wind tunnel at very low 
stagnation temperatures. 

The following section describes the calcula- 
tion method employing both eigenvalue formulations. 
Results are presented In the third section and the 
predictions of both methods are compared with mea- 
surements. The paper ends with a summary of the 
ptore Important conclusions. 

2.0 Description of the Computational Methods 

The compressible stability equations and their 
boundary conditions are well known and are given 
In several references, see for example Ref. 4. 
With the parallel flow approximation, they can be 
written In the following dimensionless form: 

Continuity: 

* dv - - dp - 

1(aU ♦ &W - w)p ♦ + 1(aU + BW) ) ♦ jjj V - 0 

( 1 ) 


♦Visiting research student from Munchen University, Germany 
♦♦Research engineer 

tVIsItlng Research student from ENSAE, Toulouse, France 
f+Assoclate Professor 
ttfProfessor and Chairman 


1 



x-Homentum 


m * 


r[1(oU * Pw - u)(<*U + ew) + (a + p 57 ) v] 

. -nJ . b 2 ) Jj . jj «. ^ . s *S) 

tH‘ dy‘ dy d 


♦ (« 2 ♦ P Z )[' 57 - 2(«u ♦ Bw) ]} 


♦ f (a 2 * P 2 )[1 £7 - (au ♦ Pw)] 


dy 


Idu, d£u dfw . T 

R dT (a 2 P p T 

R 01 dy z dy z 

1 r ,du dT x d 2 u dT du . dw., 
+ R f ( dT dy * dT 2 dy T)(o dy + P dy>> 


1 du dT r/ du A 0 dw w w 2 fl 2 .% 

* R dT d7 E(o d7 3 dy’ Ua * P )v] (2) 


Energy 
/>[1(aU + Bw - u)T 


i + d 7 ; l 


-(t - 1 >[ 


dv 


dy 

2 . »2i 


* 1(aU ♦ pw)] ♦ ^ {J-* - (a‘ + p‘)T 

* 1 r^- . d£k ,dl.2 “ . 2 dk. dl dT. 

k dT dy 2 dT 2 dy k dT dy dy 

* T * X f JL1 H l (21v(a j{ ♦ P g)v 

2 2 

du .d u d w, T 

* dT <“2 + “2> T 

dy dy 


_ii» r/„ du dw., du 

o 2 + ^2 l(a dy + P dy'' dy + 

R 

p dy' 


f „ flw du w dw_ ft dii. . . 

dy ^ dy dy * dy 


(5) 

= 0, u=>v = w = T = 0 


( 6 a) 

-> w, u, v, w t T *+ 0 


( 6 b) 


z-Homentum 


p[ 1 (aU * PW - u) (aW - PU) ♦ (a 57 - P 57) V] 


d^ „d 2 U. 


2 . „2v 


2 " p 2' 
dy^ dy^ 


11 [{a ^ - 3 *- 5 ) - (a* ♦ P fc )(aw - 0U)] 

* 2 2 
1 f du dT , dw _ du x d^ . d w Q d u . T 

* R IdT d7 t“ d7 - P d? * dT <* ~2 - P ^2 > T 


The above equations can be expressed as a first- 
order system by defining the following new vari- 
ables with primes denoting differentiation with 
respect to y. 
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In terms of these new variables, Eqs. (1) to ( 6 ), 
( 3 ) for a three-dimensional compressible flow can be 
written as _ 
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The solution of Eq. (8) subject to the boundary 
conditions given by 

y * 0, z 1 ■ z 3 ■ z 5 = » 0 

( 10 ) 

y * •, Zl . Z3* Z5 , Z7 * 0 

can be obtained with two separate eigenvalue formu- 
lations as described In the following two subsec- 
tions, 2.1 and 2.2. 

2.1 Eigenvalue Formulation Based on the Direction 
of the Amplification Vector (ONERA/CERT Method) 

The basic assumption of this method, first 
proposed by Mack/ Is that on a wing with x and 
z denoting the coordinates normal to the leading 
edge and spanwlse direction, ampll f jcatlon only 
occurs In the x-dlrectlon and not In the spanwlse 
direction; that Is, 1 


Pi - 0 (11) 

This assumption, formulated for an Infinite swept 
wing, reduces the number of unknown eigenvalues In 
the spatial theory from six to five. With u and 
R given In the transition prediction problem, the 
unknown eigenvalues correspond to , <* r and p r . 
For an assumed p r , the wave number a Is then calcu- 
lated so that the amplification rate ai can be 
determined. Additional calculations are then per- 
formed for different values of p r In order to 
determine the maximum amplification rate. This 
procedure, as In the saddle-point method, Is 
repeated for each x-statlon and the n-factor In 
the e n -method Is calculated from 

" ■ - [ <“i w x < 12 > 


for different specified frequencies. Here x 0 
corresponds to the x-statlon where the stability 
calculations are Initiated. Transition Is assumed 
to occur for the frequency for which the locally 
most unstable disturbance reaches a value of n 
between 8 and 10. 


and Is described In some detail In Ref. 7. In this 
method the relationship between the two wave num- 
bers a and p Is not assumed but computed from the 
requirement that 3a/ap Is real. According to this 
requirement, the wave orientation and growth direc- 
tion of the disturbance are given by 
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where y denotes the angle that the disturbance 
makes with the x-axIs, a and p are related by 
Eq. (13) and the disturbance propagating along the 
ray Is given by the two terms on the right-hand 
side of Eq. (13). The disturbances are damped If 
the amplification rate r defined by 


r 


- P 1 
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Is > 0, neutral If r * 0, and amplified If r < 0. 
Once a and p are computed with the constraints of 
Eq. (13), the amplification rate Is obtained from 
Eq. (14); additional calculations are then made for 
different values of 3a/3p so that new values of a 
and p are calculated to determine the maximum value 
of r. Further details of the solution procedure 
are given In Ref. 8. 

3.0 Results and Discussion 

3.1 Experimental Conditions 

The experimental data used In our studies cor- 
respond to laminar flow on a 15-degree swept 
tapered wing. The chord Is 0.228m at the root and 
0.145m at the tip. The wing has a span of 0.39m 
with an AS409 cross-section and a tralllng-edge 
sweep angle of three degrees. The height of the 
hollows on the wing Is less than 0.05 mm (from 
peak to valley) for a chordwlse extent of about 2 
cm. In order to avoid the need to perform full 
three-dimensional stabl llty/transltlon calcula- 
tions, the measurements discussed In Ref. 9 and 
summarized In Ref. 10, were carried out under 
Infinite swept conditions with the wing having a 
mean sweep angle of 12 degrees at an angle of 
attack of 0.3 degrees. The computations were made 
at a mean chord of 0.186m. 


2.2 Eigenvalue Formulation Based on the Saddle- 
Point Method fCSC Method! 

The Cebecl-Stewartson-Chen (CSC) method was 
first proposed and used by Cebecl and Stewartson 


Figure 1 shows the Mach number distributions 
measured at different Reynolds numbers for a free- 
stream Mach number of 0.74 at two stagnation pres- 
sures. As can be seen, -the Mach number distribu- 
tion has "bumps" around x/c « 0.3 and 0.47 caused 
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by small hollows In the model. On the lower sur- 
face, a smooth bend around x/c = 0,20 Is visible. 
These discontinuities are common for all experi- 
ments and their effects on transition were part of 
the Investigation conducted and discussed In detail 
In Ref. 10. 

Figure 2 shows the variation of the experi- 
mental wall temperature distributions along the 
chord. Each distribution Is referenced to Its 
temperature, Ta v , measured at the first thermo- 
couple In the flow direction. The maximum relative 
difference between the wall and freestream temper- 
ature Is 5 degrees, which Indicates that the Influ- 
ence of a nonadlabatlc wall on transition must be 
considered. This Includes the relative undulation 
as well as the absolute wall temperature. 

Figure 3 shows the variation of the drag coef- 
ficient with Reynolds numbers at three stagnation 
pressures and several stagnation temperatures. 
Since the drag coefficient Increases significantly 
for a turbulent boundary layer. It can be used as 
an Indicator of transition. The evolution of 
was measured In the cryogenic, transonic wind tun- 
nel, T2, at ONERA/CERT for different stagnation 
pressures (1.7 to 2.5 bars) and different stagna- 
tion temperatures (between 165 and 109K). The val- 
ues of C(j were obtained from wake measurements. 

In the present study we consider two runs cor- 
responding to 42 and 79, with Run 42 having a 
stagnation temperature of 145K and a stagnation 
pressure of 2.0 bar. The Influence of higher temp- 
erature at a higher pressure Is examined with run 
79, which has a stagnation temperature of 165K and 
a stagnation pressure of 2.5 bar. In both cases 
we calculate only the upper surface of the wing. 



Fig. 2. Difference between the wall and freestream 
(b) 275 bar. 



Fig. 3. Variation of the drag coefficient with 
stagnation temperature T^ and pressure p^ as a 
function of Reynolds number R c . 


Additional studies are In progress for run 60, 
which has a stagnation temperature of 7 34k, and a 
stagnation pressure of 2.0 bar. These studies will 
be reported separately. 

Figure 4 shows the experimental transition 
locations for runs 42 and 60. The location of 
transition was determined from the change of the 
wall temperature measured by thermocouples, along 
the chord, resulting from the different heat 
fluxes for laminar and turbulent flow. For run 
79, the location of transition was assumed to be 
the same as that of run 42 because of the similar 
drag coefficient and Reynolds number (Fig. 3). 
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A detailed discussion of the calculations for 
the experimental work described In the previous 
subsection Is presented In Ref. 10. In this 
subsection we present results for runs 42 and 79 
for Mach number distributions containing "bumps* 
around x/c * 0.3 and 0.47 caused by small 
hollows In the model and compare them with the 
calculations employing the CSULB method in subsec- 
tion 3.3. Studies are In progress for Mach number 
distributions without bumps and will be reported 
separately. 

The boundary-layer and stability calculations 
for run 42 were performed for an adiabatic wall 
and specified wall temperature distributions at a 
Reynolds number of 12.8 x 10 6 . Figure 5 shows 
the predictions of the ONERA/CERT method for the 
adiabatic wall case together with the measured Mach 
number distribution (Fig. 5a) used In the boundary- 
layer calculations for a stagnation temperature of 
145K and a stagnation pressure of 2.0 bar. Accord- 
ing to Fig. 4, the location of transition corres- 
ponds to x/c * 0.47. 

The calculated n-factors shown In Fig. 5b were 
obtained for seven frequencies which can be 
classified In three groups: (1) the high ones from 



Fig. 5. Distributions of (a) Mach number, (b) 
computed n factors, and (c) the most unstable 
direction of the disturbances for the adia- 
batic wall temperature distribution In Run 42 for 
the upper surface. 
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15 to 30 kHz, (2) the range from 9 to 15 kHz, and 
(3) the low ones from 2 to 9 kHz. Beginning with 
high frequencies, the calculations Indicate rather 
strong undulations of n and of the wave directions 

(Fig. 5c) along the chord. The undulations 
or n Increase with frequency and show a strong 
dependence on the pressure gradient. Excluding 
the region of high rise of H e (x/c * 0 to 0.06, 
the following mechanism can be observed. A rela- 
tive low pressure gradient has a stabilizing 
effect, whereas a higher, even positive pressure 
gradient Is destabilizing. This effect becomes 
more and more distinct with Increasing x/c. Look- 
ing, for example, at the n-curve with frequency 
30770 Hz, we can see that the relative low pressure 
gradient between x/c * 0.06 and 0.2 leads to an 
almost constant n. The Increase of dp/dx between 
0.2 and 0.3 leads to a strong rise of n. The 
steeper Mach number distribution from x/c a 0.3 
to 0.41 leads to a restabll Izatlon. The process 
Is repeated as the pressure gradient changes again 
significantly at x/c = 0.41 and 0.47. 

Figure 6 shows the variation of the computed n 
values obtained for seven frequencies at the mea- 
sured transition location of x/c = 0.47. The max- 
imum value of n Is around 4.8 for a frequency of 
approximately 25 kHz, The computed value of n for 
transition Is much lower than the value of n for 
this cryogenic wind tunnel T2 of CERT which has a 
freestream turbulence level of about 0 . 2 % and a 
transition value of n between 7 and 8 based on 
experiments performed at ambient flow conditions. 



Fig. 6. Variation of n as a function of frequency 
at the transition point x/c » 0.47. 


To investigate the Influence of wall temper- 
ature on the stability calculations, the following 
studies were conducted for specified wall temper- 
ature distributions. Figure 7 shows that the 
experimental wall temperature Is higher than the 
adiabatic wall temperature and varies more along 
the chord. Since a higher wall temperature makes 
the boundary layer more unstable. It Is plausible 
to assume that the stability calculations will 
yield higher values of n than those corresponding 
to adiabatic wall temperatures »_ 

Figure 8 shows the computed results for the 
same two high and low frequencies studies previ- 
ously. The characteristic shape of the curves has 



Fig. 7. Comparison of measured and adiabatic wall 
temperature distributions along x/c. 




Fig. 8. Distributions of (a) computed n factors 
and (b) the most unstable direction ^ of the 
disturbances for the experimental wall temperature 
distribution in run 42. 
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Fig. 9. Comparison of n-factor distributions for adiabatic and experimental wall temperature distribu- 
tions In run 42, (a) high frequency case, ( l b) low frequency case. 




Fig. 10. Comparison of distributions for adiabatic and experimental wall temperature 

distributions In run 42, (a) high frequency case, (b) low frequency case. 


not changed and the previously discussed mechanisms 
for the adiabatic wall temperature calculations are 
still valid. However, In direct comparison with 
the adiabatic case, It can be observed that the 
results with high frequency (Fig. 9a) show a 
remarkably higher n distribution for x/c greater 
than 0.17 for the case of the experimental temper- 
ature distribution while the n curves for the lower 
frequencies are almost similar (see Fig. 9b). The 
t|»m distributions for either case (see Fig. 10) 
snow that the most unstable directions are Inde- 
pendent of the wall temperature. 

Figures 11 to 13 show the results for run 79 
at a chord Reynolds number of 13.4 x IQ 6 . This 
case essentially has the same drag coefficient and 
Reynolds number as run 42. For this reason we 
assume the location of transition to be at x/c « 
0.47. 

Figure 11 shows the calculated n-factor dis- 
tribution for an adiabatic wall temperature dis- 
tribution and Fig. 12 shows the comparison of the 
n-factor distributions for adiabatic and experi- 
mental wall temperature distributions. As shown 
In Fig. 13, the temperature differences between 
adiabatic and measured wall temperatures for this 
run Is weaker than for run 42. As a result, the 
distribution of n factors for both cases do not 
differ much from each other although, as In run 



Fig. 11. Distribution of computed n factors for 
adiabatic wall conditions In run 79. 


42, the measured wall temperatures lead to higher 
values of n than those obtained with adiabatic wall 
conditions. It can also be seen that, while the 
value of n at transition location, x/c » 0.47, Is 
slightly over 7 for calculations performed with a 
measured wall temperature distribution, It Is 
around 6.3 for adiabatic wall conditions. 
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Fig. 12. Comparison of computed n factors t or adi- 
abatic and measured wall temperature distributions 
In run 79. 



We note from the above re sults that near the 
leading edge (x/c < 0.10), all the unstable fre- 
quencies are of the crossflow type. Further down- 
stream, the most unstable direction of the high 
frequencies decreases, leading to a m ore or less 
streamwlse Instability. It Is Interesting to note 
that this streamwlse Instability Is very sensitive 
to small Mach number variations (hollows) as well 
as to the wall temperature. By contrast, the most 
unstable direction of the low frequencies remains 
close to 90*: they correspond to a crossflow Insta- 
bility which does not "see* the hollows (the evolu- 
tion of the n-factor Is monotonic). This Instabil- 
ity Is also not sensitive to the wall temperature 
variations. This Is due to the fact that It Is an 
"Inflectional" Instability governed by an Inflec- 
tion point located near the outer edge of the 
boundary layer. 

3.3 Calculations with the CSC Method 

A detailed description of the stability calcu- 
lations for runs 42 and 79 are given In the previ- 
ous section for the upper surface of the AS 409 


wing. Similar detailed calculations have also been 
performed for the same runs. Including run 60, 
using the CSC method described In subsection 2.2 
and will be reported separately. In this subsec- 
tion we shall present a summary of the predictions 
of the saddle-point method of Cebecl and Stewartson 
for the same runs, 42 and 79, by showing the dis- 
tribution of n factors at several frequencies. We 
shall also present and discuss the procedure of 
determining the frequencies used In the calculation 
of amplification rates, which Is different than any 
other method which employs a combination of linear 
stability theory and e n -method to predict tran- 
sition. In fact, the studies conducted with this 
method for Incompressible flows on wings and bodies 
of revolution and recent studies In compressible 
flows over modern transport and military wings show 
that the calculation of the critical frequencies 
Is the most Important aspect of the transition pre- 
diction procedure using stability theory. The 
critical frequencies originate In a very narrow 
regions and require care and patience to compute 
their magnitudes and locations. 

The frequencies needed In transition calcula- 
tions are computed from zarfs recommended by Cebecl 
and Stewartson. 7 They essentially correspond to 
neutral stability curves In three-dimensional flows 
and have the following properties, 

«, = e 1 - o. • real (15) 

In the saddle-point method, for given velocity 
profiles obtained with the Infinite-swept boundary- 
layer method of Cebecl, the stability calculations 
begin on the zarf where, with R known and , p^ 
zero, the eigenvalue problem consists of calculat- 
ing a r , p r and « with the requirement given by Eq. 
(13). With the eigenvalues and disturbance angle 
Y of the zarf known at a specified x/c-locatlon 
and with dimensional frequency specified, the calc- 
ulations at the following x/c-statlon are performed 
to obtain a and p again with the requirement that 
da/dp Is real. This eigenvalue procedure Is then 
repeated for different values of do/ap or y to find 
the value of y for which r Is maximum at each x/c- 
statlon. This process Is repeated for each x/c- 
statlon, and n Is calculated by evaluating the 
Integral 

x 

n • - I rdx (16) 


Figure 14 shows the zarfs for run 42 for an 
adiabatic wall and measured wall temperature dis- 
tributions. As can be seen, the frequencies orig- 
inate at nearly the same location (on a vertical 
line) and vary drastically one from another. Their 
calculation requires care and patience. A paper In 
preparation will discuss our procedure for generat- 
ing them. 

Figure 15 shows the computed n factors for the 
zarfs of Fig. 14. The results show that for adi- 
abatic wall conditions, the maximum value of n for 
experimental transition location Is around 6.3; the 
corresponding value for the measured wall temper- 
ature Is 7.3. What Is more Important, however, Is 
the fact that In the latter case, the computed 
transition location agrees very well with the 
observed location considering that the n-value for 
transition In this tunnel lies between 7 and 8. 
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(a) (b) 

Fig. 14. Zarfs for run 42 for (a) adiabatic and (b) measured wall temperature distributions. 



Computed n-factors for run 42 for (a) adiabatic and (b) measured wall temperature distributions. 


Fig. 15. 

Figure 16 shows the zarfs for run 79 with wall 
conditions corresponding to adiabatic temperature 
and Fig. 17 shows a comparison between the zarfs 
obtained under adiabatic and measured wall temper- 
ature conditions. Again the steep rise In frequen- 
cies at almost one x/c-locatlon Is noted. The zarf 
calculations were performed for a very fine x/c- 
grld, since most of the frequencies .Start their 
amplification In an Interval of l/1000 th of chord 
around 0.015. Figure 16 also shows zarfs away from 
the leading edge. These zarfs have low values of 
0 r around 10“% occur In an almost zero pressure 
gradient region and do not lead to amplification 
rates that grow significantly. 

Figure 18 shows the computed n factors for run 
79 with stability calculations performed for zarfs 


In Fig. 17 under adiabatic and measured wall temp- 
erature conditions. For this flow, the n-value Is 
much higher than those In run 42. For an adia- 
batic wall, It reaches a maximum value of around 8 
and a value of around 9 for the measured wall temp- 
erature case. If we take the n value to be 7.5, a 
mld-n value of the expected n-value range for this 
wind tunnel, then transition occurs at x/c = 0.46 
for adiabatic wall conditions and x/c » 0.47 for 
measured wall temperature conditions. This com- 
pares well with Inferred transition location of 
x/c . 0.47. 

Figures 19 and 20 show a comparison between the 
calculated n-factors obtained with both methods, 
with results of CSC corresponding only to the 
disturbance that leads to transition. As can be 
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Fig. 16* Zarfs for run 79 for adiabatic wall Fig. 17. Comparison between zarfs of run 79 for 
temperature distribution. (p adiabatic and (|) measured wall temperature 

distributions. 



Fig. 18. Computed n-factors for run 79 for (a) adiabatic and (b) measured wall temperature distributions. 



fig. 19. Comparison between the predictions of CSC (denoted by circles ©) and ONERA/CERT methods for (a) 
adiabatic, and (b) mea sured wa ll temperature distribution run 42. r -._. 


seen, for x/c around 0.45, the calculations with 
the CSC method Indicate higher values of the 
n-factor. The calculations with the ONERA/CERT 
method correlate the data with n between 6.5 and 7 
while those of CSC with n between 7 and 8. 


4.0 Concluding Remarks 


Based on the studies reported here and In Refs. 
9 and 10, the following two confronts can be made. 
First, rather crude assumptions have been made In 
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Fig. 20. Comparisons between the predictions of CSC (denoted by circles©) and ONERA/CERT methods for (a) 
ablabatlc, and (b) measured wall temperature distributions In run 79. 


the computations (Infinite swept wing assumption 
with a mean sweep angle and a mean chord). Second, 
It Is possible that cryogenic conditions (tempera- 
ture fluctuations) degrade the flow quality. This 
could explain the reason why the ONERA/CERT method 
calculates n-factors somewhat lower than those com- 
puted for ambient temperature with the same stabil- 
ity code. For stagnation temperatures lower than 
those considered In the paper, Ice crystals appear 
on the model and trigger transition. 
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